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I. INTRODUCTION 



A metric of particular interest — the minimal member of the (nondenumerable) class of monotone metrics |l]J 
that can be attached to the (n 2 — l)-dimensional convex sets of n-level quantum systems is the Bures metric 
(cf. (|]). Dittmann Q established that "the connection form (gauge field) related to the generalization of the Berry 
phase to mixed states proposed by Uhlmann satisfies the source-free Yang-Mills equation *D * Duj, where the Hodge 
operator is taken with respect to the Bures metric on the space of finite-dimensional nondegenerate density matrices" 
(cf. [|o[ p. 207]). Let us also note that Bilge et al amongst others ||-|l5| (cf. @ sec. 3.f.3], 0JT1), have 

studied the properties of Yang-Mills fields in eight dimensions (in contrast to the original, principally studied case of 
four |i"s|] ). Since the three-level (n = 3) quantum systems do, in fact, form an eig/it-dimensional convex set [ p3f — the 
n x n density matrices, in general, forming (n 2 — l)-dimensional convex sets — it is of obvious interest to attempt 
to apply the analyses of Bilge et al in light of these findings of Dittmann concerning the relation between Yang-Mills 
fields and the Bures metric M. Such an effort constitutes, by and large, the substance of this communication. In 
the presentation of our findings, we will conform to the notation and conventions employed in |Llj] , as detailed in the 
appendix to that paper, although we should indicate that in a related paper po[| , Bilge replaces the notation used in 
p| for the inner product (F, F) of an S'0(n)-valued curvature form F = (F a i,) with < F,F >. We are compelled, 
however, to point out to the reader that in our computations we have not used the F a bS themselves, but rather as 
"proxies" for them, the R a b's, describing the curvature of the Riemannian manifold (cf. pll eqs. (4), (5)] p2|, eq. 
(8)]). 

Dittmann had noted || "that in affine coordinates (e. g. using the Pauli matrices for n > 2) the [Bures] metric 
becomes very complicated for n > 2 and no good parameterization seems to be available for general n" . In the first 
part of this two-part paper p4| , we attempted to fill this lacuna, by deriving the 8x8 Bures metric tensor making 
use of a recently-developed Euler-angle parameterization of the 3x3 density matrices |2^j2^]. Six Euler angles 
(denoted a,~/,a,/3,b,6), used in parameterizing an element of SU(3), are employed. Additionally, two (independent) 
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eigenvalues (Ai, A2) of the 3x3 density matrices enter into this "Schur-Schatten"-type parameterization of the 3x3 
density matrices (^,^8| . (Of course, by the requirement of unit trace, A3 = 1 — Ai — A2.) Several of the tensor elements 
were then, in fact, found to be identically zero. (The 8x8 matrix tensor decomposes into a 6 x 6 block and a 2 x 2 
one, in correspondence to the six Euler angles and the two independent eigenvalues.) However, we were not able to 
report concise symbolic expressions in [p4f for all the tensor elements. 



II. ANALYSES 

Subsequent to the issuance of Part I ]2~i|], we have further pursued the analytical matters there, managing to derive 
symbolic expressions in relatively concise form for all the entries of the 8x8 Bures metric tensor (gij), as well as its 
inverse (<? y ). Several of the entries, however, particularly two in the inverse, remain rather relatively cumbersome 
in nature, so far resisting further simplification. By way of illustration, making use of the LeafCount command in 
MATHEMATICA to measure the complexity of an expression, the largest leaf count (that is, the "total number of 
indivisible subexpressions") for any single one of the firm's is 352, while the largest for any of the g^'s is 233, aside 
from the two relatively complicated ones, having scores of 921 and 1,744. We note here that all these entries prove 
to be independent of the Euler angle a. Additionally, introducing the transformation 7 = r — a, all the entries also 
become independent of the Euler angle a. (The entry with the largest leaf count is g TT , and the second largest, 
g Ta .) Therefore, the two tensors are fully parameterizable with six, rather than eight variables. This reduction in 
dimensionality is exploited in the numerical integrations reported below. 

We noted that the two expressions (both symmetric in Ai and A2) 



A = 3 - 7Ai + 4A 2 + 7(-l + Ai)A 2 + 4A 



(1) 



and 



B = A\\ + Af (-9 + 5A 2 ) + Ai(-1 + A 2 )(-7 + 5A 2 ) + (-1 + A 2 )(2 + A 2 (-5 + 4A 2 )) 



(2) 



frequently occurred in the tensor elements. These two expressions are displayed in Figs. 1 and 2, making use of the 
transformation employed in pa]: Ai = cos 2 £1, A 2 = sin 2 £1 cos 2 £ 2 . 




FIG. 1. The subexpression A — that is (|l|), after conversion to spherical coordinates — appearing in several tensor elements 
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FIG. 2. The subexpression B — that is ([]), after conversion to spherical coordinates — appearing in several tensor elements 
For example, 

B + A(Ai - A 2 )cos26 



7 /3/3 _ 



2(-l + Ai)(-l + A 2 ) ' 

(B + A(\i - A 2 ) cos 26) esc 2 9 
-l + 2Ai + A 2 ) 2 (-l + A 1 +2A 2 ) 2 ' 



(3) 



(4) 



= gM esc 2 (3 sec 2 j3 

y 4 ^ ' 

In the two 6x6 (Euler angle) blocks of the Bures metric tensor and its inverse, the identically zero elements are 

9rb, 9r8, gab, 9a8, 9bd , 9^ , 9°^ an d Q 0,6 . 

The volume element of the Bures metric, in particular for the case n = 3, 



| ff |dadadrd/?d6d6»dAidA 2 , (6) 
where 

r— sin 26 sin 2/3 sin 2 sin 29 (Ai - A 2 ) 2 (Ai - A 3 ) 2 (A 2 - A 3 ) 2 
Vl9l ~ STO^ (A 1 + A 2 )(A 1 +A 3 )(A 2 + A 3 ) ' () 

has been determined from certain general considerations J2^,^]. We were then able, using this formula as a test, to 
numerically validate our calculations of the gij's to a very high degree of accuracy. Also, we were able to reproduce 
- using our intermediate computation of the Ricci tensor — the results of Dittmann (3^] (cf. |33|| ) for the scalar 
curvature of the Bures metric, as applied to the n = 3 case, 

2(28e 3 - 49e 2 - 9) 
s.c. = , (8) 

e 3 — e 2 

where use is made of the elementary symmetric polynomials, 

e 3 = [I A J' e 2 = Ai - A 2 + A 2 — AiA 2 - A 2 . (9) 

3=1,3 

While the convex set of n x n density matrices endowed with the Bures metric is of constant curvature for n = 2, it 



is not even a locally symmetric space for n > 2 31 (cf. |34|]). Dittmann |32| ] argues that the scalar curvature achieves 
its minimum for the fully mixed states, while others |35fc using another monotone metric (the Kubo-Mori one), argue 
that the scalar curvature is maximum for the fully mixed states. This disparity is clearly a matter of convention, 
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depending upon what sign is designated for the Ricci tensor. (Somewhat relatedly, Grasselli and Streater have shown 
that "in finite dimensions, the only monotone metrics on the space of invertible density matrices for which the (+1) 
and (-1) affine connections are mutually dual are constant multiples of the Kubo-Mori metric" |36|.) 

Additionally, it was interesting to observe that numerical evidence we adduced indicated that a certain necessary 
and sufficient condition for Riemannian connections to satisfy the Yang-Mills equation | 



(V 9 x Ric)(Y, Z) = {V a Y Ric){X, Z), (10) 

was not fulfilled in our particular situation. Here X, Y, Z denote arbitrary smooth vector fields on the manifold in 
question — that is, the eight-dimensional convex set of 3 x 3 density matrices — which we will denote by Ms, and 
V^- the covariant derivative with respect to X. 

We utilized our determination of the g^'s and g^ 's in performing the same type of calculations that have been 
shown in [fi"l| to provide upper bounds on Pontryagin numbers for SO(n) Yang-Mills fields. These bounds are 

/ (F,F) 2 >kf Pl (E) 2 , (11) 

J M 8 J M s 

where 

8 8 

(F,F) = 2^]T(F Qb ,F Qb ), (12) 

a=l b>a 

and 

(F 2 ,F 2 )>k> f P2 {E). (13) 

Here, the field strengths F^'s denote the skew-symmetric coordinate components of F, an 5 l O(n)-valued curvature 
two- form on the Euclidean space R 8 . For the zj-entry of an 8 x 8 matrix F a b was computed by us as Rij a b — that 
is, an element of the fully covariant (0,4)-Riemann curvature tensor 40, eq. (6.6.1)] (cf. jnl eqs. (4), (5)]). Also, k 



and k' are certain constants, while Pi(E) and P2(E) are the first and second Pontryagin classes for the bundle E in 
question. The principal U (n)-bundle E (n = 3) utilized by Dittmann || is Gl(n, C) — > Gl(n, C)/U(n), where Gl(n, C) 
denotes the general linear group over the complex numbers C, that is the nonsingular n x n matrices with complex 
entries J4l| , ex. 5.1]. (The Pontryagin number for a /cmr-dimensional compact manifold M4 is equal to the integral 
over M4 of a representative of the first, and unique Pontryagin class of the bundle, that is the Chern number j M 72 
of the complexified bundle pS^] .) 

The first of the two bounds ( |Tl| ) was reported in jllj] and the second ( |l3| ) in Jl5| . The first is achieved by strongly 
[anti-] self-dual Yang-Mills fields, while the second is too restrictive in this regard 111]]. An additional curvature 
invariant is (trF 2 ,trF 2 ), so a "generic action density" can be taken to be of the form jO, eq. (f7)] 

p(F, F) 2 + q(F 2 ,F 2 ) + r(trF 2 , trF 2 ). (14) 

("The action density should be written in terms of the local curvature 2-form matrix in a way independent of the 
local trivialization of the bundle. Hence it should involve invariant polynomials of the local curvature matrix. We 
want to express the action as an inner product in the space of fc-forms, which gives a quantity independent of the 
local coordinates" fill] .) 



A. Monte-Carlo numerical integrations 

To evaluate the integrals over Mg of the three curvature invariants occurring in the generic action (pT[), we pursued 
a Monte-Carlo strategy Jil^-fi^] in which we randomly selected points in the six-dimensional hyperrectangular domain 
determined by the six active parameters (four Euler angles — r, /3, b, 9 — and two independent eigenvalues reexpressed 
using spherical coordinates, that is, Ai = cos 2 £1, A2 = sin 2 Ci cos 2 (2 |§^])- This six-dimensional domain has one side 
of length 7r (corresponding to the range of r), three sides (for (3, b and 9) of length tt/2, one side (for (1) of length 
cos -1 7r/3 ~ .955317 and one (for £2) of length 7r/4. pq ]. The two remaining (seventh and eighth) sides — associated 
with the Euler angles a and a, which we have noted are absent from the simplified transformed expressions for the 
gij's ]^4| — are of length n. At each of the randomly selected points, we computed the products of the volume 
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element (§) with the quantities (F, F) 2 , {F 2 ,F 2 ), (trF 2 ,trF 2 ) — as well as (F 3 ,F 3 ) 2 / 3 and (T 4 ,^ 4 ) 1 / 2 . (In these 
MATHEMATICA computations, we utilized the relation det(I + tF) = X]fc=o <Jfc ^ fe ' wnere (J k are invariants of the 
local trivialization [jll], eq. (12)], Q.) We multiplied the averages of each of the five products over the randomly 

7 11 

chosen points by the full eight-dimensional Euclidean hyperrectangular volume, which is ^ cos 775 ~ 90.1668. For 
21,238 randomly generated points, our set of five results (exactly as computed, without any "rounding") is presented 
in the first line of Table I. 



field 






J Ms (trF 2 ,trF 2 ) 


j Ms (F\ F r- /s 


j Mx (F\ F r /2 


Bures metric 


.00174878 


9.91872 • 10~ 14 


.00174878 


7.90173 • 10" 22 


5.07518 ■ HT™ 


anti-self-dual part 


2.0692 


1.5519 


3.62111 


1.41 


1.09736 


self-dual part 


18.6228 


7.75951 


26.3823 


5.15961 


3.29208 



TABLE I. Curvature invariants approximated by evaluating the integrands at 21,238 points randomly selected in the 
six-dimensional hyperrectangle of active parameters 
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We also, following Baulieu and Shatashvili 0, eqs. (3.4)-(3.6)] (cf. @ p. 207], @, eq. (2.9)]), decomposed 
the curvature two-form for the 8x8 Bures metric tensor into two S'pm(7)-irreducible components, according to 
28 = 21 ® 7, which can be called self-dual and anti-self-dual respectively. We then evaluated the same curvature 
norms. The corresponding Monte-Carlo results for the same set of 21,238 randomly selected points are presented in 
the second and third lines of Table I. It clearly appears, then, that the Yang-Mills field for the Bures metric on the 
three-level quantum systems is neither "self-dual" nor "anti-self-dual" (cf. ]49|-|53|] ) . 

Since all the five integrands used in Table I are nonnegative in nature, it appears that for the Yang-Mills/Bures 
field (corresponding to the first line of the table), (F 3 ,F 3 ) and (F 4 ,F 4 ), at least, are themselves zero. Also, we see 
that in in this same first line (but not the second and third), the integral of (F, F) 2 equals that of (trF 2 , trF 2 ). Here 
eq- (21)] 
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trF 2 = -2j2J2 F « b - ( 15 ) 



a—l b>a 



(For the octonionic 50(7) instanton two-form F, the four-form TrF A F is neither self-dual nor anti-self-dual f34j.) 
The 8x8 skew-symmetric matrices F a b must have four pairs of imaginary eigenvalues of opposite sign. ("Strong self- 
duality or strong anti-self-duality can be characterized by requiring the equality of the absolute values of the [eight] 
eigenvalues" [jll).) For the eight-dimensional Yang-Mills field F over the three-level quantum systems, our analyses 
indicate that one of these four pairs is always (0,0) — clearly indicative of some form of degeneracy. (Riemannian 
manifolds the skew-symmetric curvature operators of which have constant eigenvalues have been the subject of study 
p5[ , including the case of eig/it-dimensional manifolds [Q.) This ensures that the determinant (the product of the 
eigenvalues), as well as the (F 4 b , F^ b )'s themselves, and thus (F 4 , F 4 ), are zero. The apparent zero nature of (F 2 , F 2 ) 
and (F 3 , F 3 ) can not be directly explained in this manner, that is, by the zero nature of one pair of eigenvalues. (Let 



us point out here the possible relevance for our study, yielding a "small action", of |57|.) 



B. Numerical integrations using regular lattices 

As somewhat of a "cross-check" on the Monte-Carlo analyses reported above in Table I (motivated, in part, by 
some apparent instabilities in early, preliminary analyses), we also undertook numerical integrations by evaluating the 
various curvature quantities not at randomly selected points in the six-dimensional hyperrectangle, but at the nodes 
of regular lattices superimposed on this hyperrectangle. For a 2 6 = 64-point lattice (dividing each parameter range 
into two equal segments and taking the midpoints of each segment as coordinates), we found 



field 






(trF 2 , trF 2 ) 


^^3^3)2/3 




Bures metric 
anti-self-dual part 
self-dual part 


.00109779 
1.77679 
15.9911 


1.399 ■ 10~ 2b 
1.3326 
6.66298 


.00109779 
3.10939 
22.6541 


4.297 • 10~' io 
1.21074 
4.43048 


7.231 ■ 10~ 44 
.942287 
2.82686 


SD part minus ASD part 


28.4287 


21.3215 


49.7502 


19.3719 


15.0766 



TABLE II. Curvature invariants approximated by evaluating the integrands at the nodes of a 64-point lattice imposed upon 
the six-dimensional hyperrectangle of active parameters 
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A fourth line is included in Table II for the difference jl7[ eq. (3.4)] 

= F + - F~, (16) 
of the self-dual (F + ) and anti-self-dual (F~) fields, so 

F ± = i(F±tF). (17) 

For a (finer) 3 6 = 729-point lattice, we obtained 



field 






J Ms (trF 2 ,trF 2 ) 


/mr(F 3 iF 3 )2 /3 


f (F 4 .F 4 ) 1/2 


Bures metric 
anti-self-dual part 
self-dual part 


.00029131 
2.11051 
18.9946 


1.22054 • 10~ ab 
1.58288 
7.9144 


.00029131 
3.69339 
26.909 


1.15258 ■ 10 _3 ° 
1.43814 
5.2626 


1.17924 ■ 10~ 42 
1.11927 
3.3578 


TABLE III. Curvature invariants approximated by evaluating the integrands at the nodes of a 729-point lattice imposed 
upon the six-dimensional hyperrectangle of active parameters 

Additionally, for a 4 6 = 4096-point lattice, we computed 


field 






Jm>F 2 , trF 2 ) 


^^3^3)2/3 




Bures metric 
anti-self-dual part 
self-dual part 


.000824594 
2.08472 
18.7625 


5.54578 ■ 10"™ 
1.56354 
7.81771 


.000824594 
3.64827 
26.5802 


3.95448 • 10~ JO 
1.42057 
5.19831 


9.09798 ■ 10~ 4J 
1.10559 
3.31677 



TABLE IV. Curvature invariants approximated by evaluating the integrands at the nodes of a 4096-point lattice imposed 
upon the six-dimensional hyperrectangle of active parameters 
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We see that the results here, in their overall aspects, are supportive of our much more extensive (and thus we are 
inclined to believe more accurate) Monte-Carlo analysis in Table I — to which we are in the process of adding still 
more randomly selected points. (It might be noted that none of our tabulated results — for the anti-self-dual and 
self-dual components — when divided by 7r" , for n = 4 or other small integers, gives indications of being integral in 



nature 58, eq. (7)]. Also, there might seem to be a question of whether our tabulated results should be "corrected" 
by a factor of 3! = 6, that is the number of ways of permuting the rows and columns of a 3 x 3 density matrix, but 
doing so would lead to density matrices that are simply permutations of one another being represented by distinct 
points in the manifold ^5|, |3(], sec. III.B].) 

C. Computation of Yang-Mills action functionals 

For a G-bundle E— > X, with connections A, the Yang-Mills functional can be expressed in the eight-dimensional 
case as pi, eq. (2.4)] 

S YM [A] = ^ d 8 *^ V" = ^ll^l| 2 , (18) 

where \ represents the coupling constant of the gauge fields. For the Yang-Mills gauge field defining the Bures 
metric B], based on 14,151 randomly chosen points in the six-dimensional hyperrectangle of active parameters, we 
obtained a value of .0145485 for ||-F|| 2 . For ||-F~|| , the comparable value was much larger, that is 5.33255 • 10 6 , 
and almost identical but slightly larger, that is 5.33268 • 10 6 , for ||.F + || . "It can be shown that any self-dual 2- 
form defined by the above criteria [for self-duality in dimensions greater than 4] satisfies the Yang-Mills equations. 
However, the corresponding Yang-Mills action need not be extremal" [[h| p. 4804]. ". . . it is no longer true that a 
configuration satisfying our duality equation [in arbitrary even dimensions] automatically corresponds to a minimum 
of the [Yang-Mills] action" pi sec. 2.4]. 
Also Gao and Tian assert |48|, eq. (2.5)] 

Sym[A}= 1 (2 / OATr(FAF) + v/B 2 + 8A||F+|| 2 ), (19) 

g z (B + vB Z + 8A) Jx 

where @ eq. (2.8)] A = 6,B = -4 and the four-form Q is a "bispinor" p§ eq. (2.7)]. They state that "the 
Yang-Mills action can be written as a non-negative term proportional to \\F + || , plus a topological invariant. Clearly, 
such an action will reach its minimal values at F + = 0" p8L p. 3] (cf. [^9]-^). 

In an exchange with Y. H. Gao on this matter, he has written: "Tentatively, the puzzle you mentioned (i. c. the 
Yang-Mills functional for the (anti) self-dual field is larger than for the full field) could stem from the problem of 
constructing a topological term using numerical methods. In the analytic world the difference between ||F|| 2 and 
4||i ? + || 2 is 'topological' in the sense that its variation with respect to the connection vanishes. So, locally a minimum 
of Il-F+H 2 is also a minimum of the original Yang-Mills action, and \\F\\ 2 should be larger than ||.F_|| 2 provided 
F and F_ belong to the same 'topological class' (namely the topological term has the same value). In numerical 
computations, however, manifestation of topological invariance for an integral might not be as easy as in the analytic 
world. It could be the case that when F is changed into F_, the value of the "topological term" is also changed. In 
that case the argument in the analytic consideration would be no longer valid. Thus a naive suggestion is that during 
comparison between ||F|| 2 and H-F-H 2 numerically, one should also monitor possible changes in the topological term, 
to rule out the possibility that F and F- are not really in the same topological class." 

So, it appears that it would be of interest in this matter to, in addition to our other computations, attempt to 
numerically approximate the term 



/ 0, A Tr(F A F), (20) 
Jx 



in (19) by itself, as well as with F replaced in it by F + as well as F_. 



D. Monopole equations with 5 , pin(7)-holonomy 

Motivated by the Seiberg-Witten equations p3], Bilge et al formulated eight-dimensional monopole equations as 
follows [P> eq. (24)]: 
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£U(¥) = 0, p+(F+) = (¥**)+ (21) 

Here (in the notation of pM) is the orthogonal projection of onto the spinor subbundle spanned by 

P + (/j)j* = 2 ... 7. We numerically solved the second set of equations ( pl| ) for the components of the spinor |16, eqs. 
(3a.57), (3f.l02), (3f.l03)] (cf. g p. 30]) 

* = i (ei + e 4 i, e 2 + e 5 i, e 3 + e 6 i, e + e 7 i, e x - e 4 i, e 2 - e 5 i, e 3 - e 6 i, e - e 7 i), (22) 

. Only for what we have termed the "anti-self-dual" component (F~) did we obtain solutions, which always came in 
one of two forms. The first was e 2 = ~, with u equal to a real constant, and the other six e's all set to zero, while 
the second solution was eo = e 3 = e§ = e-j = 0, 

wies uie4 

ei = -5— — 21 e 2 = "T^ — 2' ( 2 ^) 
e\ + ef ef + ej 

where u is the same as in the first solution. 

It remains for us to attempt to solve the first of the two sets of equations in that is the Dirac equation, 

together with the Coulomb gauge condition. 

III. DISCUSSION 

What we have considered to be the "anti-self-dual" part (F~) of the Yang-Mills field over the eight-dimensional 
convex set of three-level quantum systems satisfies the "set 6" of twenty-one equations in jL2| , 

F\2 — -F34 = 0; F\2 — = 0; F12 — -FYs = 0; F13 + F24 = 0; F13 — F57 = 0; F13 + Fes = (24) 
Fm - F 23 = 0; F14 + F 67 = 0; F 14 + F 5S = 0; F 15 + F 26 — 0; F 15 + F 37 = 0; F 15 — F 48 = 
Fie — F25 — 0; Fi6 — F38 — 0; Fie — F47 = 0; Fn + F2& — 0; Fu — F35 = 0; Fn + F^e = 

Fig — F27 = 0; Fig + F 36 = 0; Fig + F 45 = 

that is, to the negative eigenvalue, -3, of a fourth rank tensor invariant under SO(7). (A skew-symmetric matrix 
satisfying the set ( p4| ) must have all its [imaginary] eigenvalues equal in absolute value pd|.) 

The "self-dual" part (F + ) of the Yang-Mills/Bures field satisfies the "set a" of seven equations |l2| , 

F12 + F34 + F56 + F78 = 0; F13 — F24 + F57 — = 0; F14 + F23 — — F58 = 0; F15 — F26 — F37 + F48 = 

(25) 

Fi6 + F25 + F 38 + F47 = 0; F17 — F 2 s + F35 — F46 = 0; Fig + F27 — F 36 — F45 = 

that is to the positive eigenvalue 1. 

"The solutions of set a and set b can be viewed as analogues of self-dual 2-forms in four dimensions from different 
aspects. The strongly self-dual 2-forms, hence the solutions of set b saturate various topological lower bounds . . . , but 
they form an overdetermined system . . . we show that the solutions of set b for an TV-dimensional gauge group, depend 
exactly on N arbitrary constants, provided that the system is consistent. Thus the set b lacks the rich structure of 
the self-duality equations in four dimensions. On the other hand, the solutions of set a do not saturate the topological 
bounds . . . but these equations form an elliptic system under the Coulomb gauge condition" fll^ . 

"In fact, one often finds in the literature that [our] equations (^5|) and (E4j) are referred to as the self-duality and 
anti-self-duality equations respectively. This nomenclature suggests a symmetry between these equations which is not 
present in the octonionic case since, for example, the spaces have different dimension. In our opinion, self-duality and 
anti-self-duality correspond to which way the division algebra A acts: if on the left or on the right, and are hence 
related by a change of orientation on the manifold. Although there has been some work in the literature concerning 
equation (E4), we believe this equation is not as fundamental as (pq). This can also be seen not just in the results 
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of the present paper but also, for example in |62|| , where it is shown that supersymmetry singles out equation (Ul 
p3j. (In our examination of various papers, it even appears that what one author calls the "self-dual" part, another 
author may call the "anti-self-dual" part — for example, | ]To| , ^8| .) 

Of obvious research interest would be companion studies to that (eight-dimensional one) here of the Bures metric 
on various /owr-dimensional convex sets of density matrices — such as those examined in 0,|64| . In line with this we 



have sought to calculate the Yang-Mills functionals, as in sec. II C, for the /cwr-dimensional convex set of 3 x 3 density 
matrices studied in fl6^1 endowed with the Bures metric 

x — iy \ 

2~2v (26) 
v - z J 

The Yang-Mills functionals for the original and self-dual and anti-self-dual fields turn out to be infinite in this case, 
though. 

In the context of D-dimensional Euclidean gravity, Acharya and O'Loughlin have defined a generalization of the 
self-dual Yang-Mills equations as duality conditions on the curvature two- form R of a Riemannian manifold 




For D = 8, solutions to these self-duality equations are provided by manifolds of Spin(7) holonomy |65-67|, which 
are necessarily Ricci-flat. The holonomy group is then 21-dimensional. Acharya and O'Loughlin replaced the field 
strengths in the self-dual Yang-Mills equation 

= -(f>LLv\pF\p (27) 

by the components R a b of the curvature two- form (as we have above). 

Using the technique of embedding the spin connection in the gauge connection (cf. |]68| , sec. 10.5]), they construct 
a self-dual gauge field directly from the self-dual metric. They let u G a b be the generators of one of SU(2), SU(3), G2 
or Spin(7). The ansatz for the gauge field is A = "fG a bt-u a b. The form index of A comes from the form index of uj 
while the Lie algebra structure of A comes from that of G. One sees then that 

F = dA + A A A (28) 



= "fGabduJab + ^ 2 G ab UJ ac A LU cb (29) 

where k is a constant that depends upon the group generated by G a b- In each case it is trivial to solve for 7 giving 
F = r yG a bRab- Duality of F follows from that of R and the symmetry of R^xp between the first pair and second pair 
of indices" We shall attempt to follow this prescription in our continuing research on the topic of this paper. 

The duality operator, 4> a bcd for D = 8, is the unique Spm(7)-invariant four-index antisymmetric tensor which is 
Hodge self-dual. Proposition 10.6.7 of |67j relates the first Pontryagin class of a manifold M with a Spin-(7) structure 
(f2, g) to the integral over M of |i?| 2 , where R is the Riemann curvature of g. Also, if (M, f2, g) is a compact Spin(7)- 
mainfold, then g has holonomy Spin(7) if and only if M is simply connected and the Betti numbers of M satisfy 
6 3 + b\ = b 2 + 2bi + 25 §7|, Thm. 10.6.8]. 

It has been known for some time that the Bures metric on the two-level quantum systems is isometric to the standard 
metric on the three-sphere |4|]69]]. "The Bures metric for a two-dimensional system corresponds to the surface of a 
unit four-ball, i. e., to the maximally symmetric three-dimensional space of positive curvature (and may be recognized 
as the spatial part of the Robertson- Walker metric in general relativity). This space is homogeneous and isotropic, 
and hence the Bures metric does not distinguish a preferred location or direction in the space of density operators. 
Indeed, as well as rotational symmetry in Bloch co-ordinates (corresponding to unitary invariance), the metric has a 
further set of symmetries generated by the infinitesimal transformations 

r -> r + e{l - r 2 ) J / 2 r (30) 

(where r is an arbitrary three- vector)" [^9], p. 128], r being the radial distance — the length of r — in the Bloch 
sphere of two-level systems (from the fully mixed state, r — 0). Petz and Siidar observe that "in the case of the 
[Bures] metric, the tangential component is independent of r" p. 2667]. 

Contrastingly, Dittmann plj ] established that the Bures metric on the three-level quantum systems is not a space of 
constant curvature nor even locally symmetric. (Neither, is the Einstein- Yang-Mills equation fulfilled with a certain 
cosmological constant M, as it is for the two- level quantum systems |7Q].) Nevertheless, we have obtained in Tables 
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I- IV, strong indications that the curvature (known to satisfy the Yang-Mills equation ||) of the Bures metric on the 
three-level quantum systems is quite flat in character. ( "quantum information manifolds are equipped with two natural 
flat connections: the mixture connection, obtained from the linear structure of trace class operators themselves, and 
the exponential connection, obtained when combinations of states are performed by adding their logarithms... the 
Bogoliubov-Kubo-Mori metric is, up to a factor, the unique monotone Riemannian metric with respect to which the 



exponential and mixture connections are dual |36|.) 

Due to the relative complexity (high leaf count) of our current formulas for the entries of the 8x8 Bures metric 
tensor and its inverse, we have had to have recourse here, by and large, to numerical calculations in order to gain 
insights into the nature of the associated Yang-Mills field. If further progress in simplifying these entries can be 
achieved, it may be possible to proceed further with more highly desirable exact symbolic computations. 
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